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> ; Abstract 

1-2 ' In a broad class of dark energy models, the universe may collapse within a finite time tc- Here 
we study a representative model of dark energy with a linear potential, V{(l)) = Vo(l + acj)). This 
■ model is the simplest doomsday model, in which the universe collapses rather quickly after it stops 
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expanding. Observational data from type la supernovae (SNe la), cosmic microwave background 
anisotropy (CMB), and large scale structure (LSS) are complementary in constraining dark energy 



■ models. Using the new SN la data (Riess sample), the CMB data from WMAP, and the LSS data 

■ from 2dF, we find that the collapse time of the universe is t^. > 42 (24) gigayears from today at 68% 
Q"^ ■ (95%) confidence. 

o . 

^ ■ PACS numbers: 98.80.Cq, ll.25.-w, 04.65. +e 
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I. INTRODUCTION 



The existence of dark energy 0, 0] is one of the most significant cosmological discoveries 
ever made. The nature of this dark energy has remained a mystery. Viable models of dark 
energy include a cosmological constant, quintessence A,i5||, or modified Friedmann equations 
0. The only available model of dark energy based on string theory describes a metastable 
cosmological constant with an exponentially large decay time . See for reviews with more 
complete lists of references. 

Although current observational data seem to favor a cosmological constant 0, 11, 1^ . 



there are two good reasons to continue dark energy studies. First, the measured dark energ 
density^ as an arbitrary function of time, X{z) = px{z)/ px{0), has large uncertainties [1( 



aensity 

3113: 



hence it does not rule out dark energy models with px{z) that do not vary greatly 

with time Second, even if px{,z) were measured to be a constant at high accuracy in 

future observations, we still need to find physically well motivated models that give an effective 
cosmological constant with an incredibly small value A ~ 10~^^ g/cur". 

One of the first (and by far the simplest) models of dark energy, which we are going to discuss 
in our paper, was proposed in |3| in order to address the cosmological constant problem. In 
this model the cosmological constant was replaced by the energy density of a slowly varying 
scalar field with the linear effective potential 

n0) = yo(l + «0). (1) 

Here we use the units Mp = (SttG)"^/^ = 1. If the slope of the potential is sufficiently small, 
c^Vq ^ 10~^^° in Planck units, the field (p practically does not change during the last 10^° years, 
its kinetic energy is very small, so its total potential energy V{(f)) acts nearly like a cosmological 
constant. The anomalous flatness of the effective potential in this scenario is a standard feature 
of most of the models of dark energy 0, 0, 1^- 

Even though the energy density of the field in this model practically does not change at the 
present time, it changed substantially during inflation. Since is a massless field, it experienced 
quantum jumps with the amplitude H/2tt during each time interval H^^ . These jumps moved 
the field in all possible directions. As a result, the universe became divided into an infinitely 
large number of exponentially large 'universes' containing all possible values of the field 0, 
and, consequently, all possible values of the effective cosmological constant A(0) = 1^(0). This 
quantity may range from —M^ to +Mp in different parts of the universe, but we can live only 
in the 'universes' with |A| < 10^^^ g/cm^. 

Indeed, if A < —10^^^ g/cm^, the universe collapses within a time much smaller than the 
present age of the universe ~ 10^" years 0, 0, 3]- On the other hand, if A ^ 10~^^ g/cm^. 



the universe at present would expand exponentially fast, energy density of matter would be 
exponentially small, and life as we know it would be impossible 0, 0| . This means that we 
can live only in those parts of the universe where the cosmological constant does not differ too 
much from its presently observed value |A| ~ 10~^^ g/cur' . This approach constituted the basis 



^ The dark energy densit y is better constrained by data than the dark energy equation of state, since it is more 
closely related to data.|l3l| 
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for many subsequent attempts to solve the cosmological constant problem using the anthropic 
principle in inflationary cosmology [1 



There are several other models of a slowly rolling scalar field which may allow one to solve 
the cosmological constant problem. For example, one may consider models with quadratic 
potential rv?^ jl — A 3. . .20.1 . with potentials Vi^e°"^ — A ji^], potentials unbounded from below 
in = 8 supergravity |22j |. or potentials depending on two different fields aVo0 + fi^ jl^. 
In all of these models the anthropic solution of the cosmological constant problem is based on 
the assumption that the effective potential continuously changes in a sufficiently broad interval 
near = 0. In all of these models the scalar fields eventually roll down to negative values of V ^ 
and the universe experiences global collapse, just as the universe with a negative cosmological 
constant. A detailed description of this process can be found in jisi ]. 

Thus, in all of the models mentioned above, the anthropic solution of the cosmological con- 
stant problem goes hand in hand with the prediction of the global collapse of the universe.'^ This 
makes it especially interesting whether there could be any observational evidence in favor of 
these models, and what kind of bounds on the lifetime of the universe one could obtain by 
cosmological observations. 

Talking about the future observations, even if the results to be obtained by JDEM and 
Planck satellites will favor the simplest cosmological constant model, they will not necessarily 
imply that the universe is going to expand forever. They will only place some model-dependent 
limits on the possible time of the doomsday. For example, in the context of the simplest linear 
model (^, such data would only imply that the collapse of the universe cannot happen earlier 
than tc ~ 40 gigayears from now, at the 95% confidence level [s^ . 

The actual data to be obtained many years from now may alter these estimates. It is 
important to understand what we can say about the possible doomsday time on the basis of the 
presently available data. There were several attempts to study this question 



but the results of these investigations appear to be strongly model-dependent. Our approach 
will be similar to that of Ref. where the constraint on the doomsday time tc > 30 Gyrs 
was obtained for models with the potential Voe'^^'^ — A. However, this constraint would change 
significantly if one considered a more general function Voe°"^—A depending on three parameters, 
a, Vq, and A, in addition to the initial value of the field 0. 

In this paper we are going to study the observational constraints on the lifetime of the uni- 
verse in the simplest but rather general class of models where the effective potential with respect 
to the slowly rolling field can be represented by the linear function ([Q) in the cosmologically 
interesting interval ^ 10~^^ g/cm^ An advantage of this model (see also [3^) is that 

it depends only on one parameter aVo- This can be made manifest by a shift — >■ — Vo/a 
which reduces the potential (1) to an even simpler expression q;Vo0. This will allow us to obtain 
more definite constraints on the doomsday tc- 

We will use the Markov Chain Monte Carlo (MCMC) technique, which will allow us to 
compare observational data with the model ([Q) with different values of a. Even for this sim- 



^ One can avoid this conclusion in the models where the cosmological constant can take an exponentially large 
number of discrete values near V — EH Ell • In the string theory based models of this type the universe 
may either collapse or spontaneously decompactify and become 10-dimensional 
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pie model, we still need an additional input from the theory. The simplest and perhaps the 
most reasonable approach is to say that since we do not know anything about the origin of 
the parameter a, we will consider all values of a being equally probable, and use fiat prior 
probability distribution = const as an input in the MCMC calculations. But various 
theoretical arguments may give different prior probability distributions p{a), which may affect 
the results. For example, one may argue that if there is some reason why the parameter aVo in 
is smaller than 10~^^°, in Planck units, then perhaps it is many orders of magnitude smaller 
than 10~^^° js^]. This is a very strong argument which suggests that p{a) may be strongly 
peaked at small a, in which case the theory will have almost exactly the same observational 
consequences as the simple cosmological constant model. On the other hand, instead of the 
model (0) one may consider models with potentials aVo</> + /9C". During inflation, each of these 
two fields will fluctuate, so that in different parts of the universe with \V\ < 10"^*^ g/cm^ one 
will have different combinations of the fields and ^, and the slope of the effective potential 
in the steepest descent direction will also take different values. Such models in the observa- 
tionally important region \V\ < 10~^® g/cm^ behave just as the model with a constant slope 
of the effective potential (Q) with respect to a combination of the fields (j) and but the prior 
probability distribution to live in a part of the universe with a given effective slope aVo will 
depend on the parameter n ji^ . 

That is why in our calculations we will examine the linear model using the simplest 
prior probability p{a) = const, as well as two other possibilities, p{a) oc a and p{a) oc 
Whereas the resulting constraints on the lower bound on the doomsday time tc appears not to 
be very sensitive to the choice of p{a), one should keep in mind that, in general, predictions of 
the fate of the universe depend not only on the choice of the potential, but also on the prior 
probability for the parameters of the model. This does not make the results of the MCMC 
calculations ill defined, it just means that the definition of a physical model of dark energy 
should include some idea about the probability of the parameters of the model, as well as the 
probability of initial conditions for the scalar field. In the absence of any a priori idea, it is 
reasonable to use a fiat prior for these parameters. 

In Sec. II of our paper we will describe our analysis technique and the data used. We present 
our results in Sec. III. In Sec. IV we will make some methodological comments, comparing the 
results of the MCMC calculations with the Fisher matrix approach. Sec. V contains a critical 
discussion of the possibility to obtain model-independent constraints on the doomsday time. 
We summarize our results and discuss their implications in Sec. VI. 



II. METHOD 

A. Evolving the scalar field 

Following Ref.js^l, we assume a fiat universe, Qtot = + = 1, which is consistent with 
current observational data 0, ■ Note that since the expansion rate of the universe today, 
Hq, depends on Qm, Vq and a, the only independent parameters in this model are {flm, a, -f^o)- 
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The usual equations of motion are 



+ 3H{a)(f) + 
SttG ' 



dV 



(2) 
(3) 



where a is the cosmic scale factor, is the matter density today, and the dots denote derivatives 
with respect to the cosmic time t. The Hubble parameter is 



H''{a) 



SttG 



a-^ 2 



(4) 



It is not straightforward to solve the usual equations of motion for a given set of (fi^, a, 
Hq), since flm depends on Hq, which in turn is only known once we have solved the equations 
for a{t) and found its derivative at a = 1. 



We rewrite the equations of motion as 

— + 3^4^) -i- + a 

(Jj I Lt I 



1 d'^a 



1 + ac 



1 da Y 
aar J 



Eliz) = 





1 n* 
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d^j 


' 2 a3 




1 /d0 







dr 



+ 1 + ac 



where 



8nG 



a 



(5) 
(6) 

(7) 
(8) 



3 ' " V 3 
Note that the equations of motion (Eqs.[S]-[7]) only depend on {fl'^, a). 

At r — s> 0, a — > 0, we take (p = d(j)/dt = 30]. Evolving the equations of motion to a = 1 



0) gives us 



E,{z = 0), 



EH 



0)' 



(9) 



B. Prior distributions of parameters 



In the linear model of Eq.([P), the expansion rate of the universe depends on the parameters 
(Pm' ^) (s6e Eq.(jll)). We assume uniform priors for and Vq. It can be shown that 
choosing uniform prior for Vq is equivalent to choosing uniform prior for the initial value of 
the scalar field after inflation, which is the standard assumption made in 0,1^. For a, we 
consider three different priors (see Ref.ji^ for justifications): 

P{a) oc 1; P{a) oc ^ P{a) oc a (10) 
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which correspond to n = 2, 3, 3/2 in Eq.(29) of Ref. 2^ . 



Using the Bayes' Theorem, we can take the prior in a into consideration by replacing the 
usual with 

(11) 



x' = x'-2 1nP(a). 



C. Data used 



The measured distance-redshift relations of type la supernovae (SNe la) provide the founda- 
tion for testing dark energy models. In a fiat Universe, the dimensionless luminosity distance 

dL[z) = cE^^(\^z)V{z), V{z)= [ dz'/E{z'), (12) 

Jo 

where 

E{z) = E,{z)/E,{z = 0) (13) 

with E^:{z) given by Eq.(|7j). 

The observational data we use in this paper are the same as in • We use the "gold" set of 
157 SNe la (the Riess sample) published in and analyze it using fiux-averaging statistics [l^ 
3^ to reduce the bias due to weak gravitational lensing by intervening matter p^. A Fortran 
code that uses fiux-averaging statistics to compute the likelihood of an arbitrary dark energy 
model (given the SN la data from |9|]) can be found at hUp://www.nhn.ou.edu/r^wang/SNcode/ 



We only use CMB and LSS data that are not sensitive to the assumptions made about dark 
energy 



inly use 

la 111. 



The only CMB data we use is the measurement of the CMB shift parameter |3^ 

R = n]l^T{zcMB) = 1.716 ± 0.062 (14) 

from CMB (WMAP, CBl, ACBAR) HliJ, where zcmb = 1089. 

The only large-scale structure information we use is the linear growth rate f{z2dF) = 0.51 ± 
0.11 measured by the 2dF galaxy redshift survey (2dFGRS) |4l[ Hi, where Z2dF — 0.15 is the 



effective redshift of this survey. Since f = f3 * bi, where /3 is the redshift distortion parameter 
measured from the ratio of the redshift-space correlation function to the real-space correlation 
function [see Eq.(17) in j4l|], and 6i is the bias factor [square root of the ratio of galaxy power 
spectrum and mass power spectrum]. Since both correlation function and power spectrum are 
statistical descriptions of galaxy survey data, they can be extracted from data without making 
specific assumptions about dark energy. Note that on the other hand, the theoretical prediction 
of the linear growth rate does depend on assumption about dark energy, as well as cosmological 
parameters; this is why we can use galaxy survey data to probe dark energy and constrain 
cosmological parameters. For a given set of cosmological parameters and an assumed dark 
energy density px{z), the linear growth rate / = (cilnD/cilna) is determined by solving the 
equation for the linear growth rate D, 

D"{t) + 2E{z)D\r) - ^fi„(l + zfD = 0, (15) 
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where primes denote d/d{Hot), and E{z) = + z)^ + (1 - ^m)px{z)/px{0)]^^^ for a flat 

universe. For the dark energy model studied in this paper, E{z) is given by Eqs.lfT^ and ((Tj). 

We use the Markov Chain Monte Carlo (MCMC) technique in the likelihood analysis (based 
on the MCMC engine of ^^), and obtain a few million samples of (i^^; ^5 -^o)- This method 
samples from the full posterior distribution of the parameters, and from these samples the 
marginalized posterior distributions of the parameters (i.e., their probability distribution func- 
tions [pdf]) can be estimated. We use the method proposed by Gelman and Rubin to test for 
convergence 0,13. This method uses a convergence indicator 

[{N-1)/N]W + B^{1 + 1/M) 
R = ^ , (16) 

where M is the number of chains (each with 2N elements) starting at well-separated points in 
parameter space, W is the mean variance of the chains, and is the variance between the 
chains. Convergence is achieved for R < 1.2. We find that R < 1.01 for our MCMC chains, 
which assures us that convergence has been achieved. 



III. RESULTS 



In models that lead to a Big Crunch, the Hubble parameter H{z) oc E{z) will decrease with 
time until E(z) = at t = ttum, when the universe stops expanding and starts to collapse. 
Fig.l of Ref . |3jj| shows the cosmic scale factor in five models (the linear potential model with 
five different parameter choices). Clearly, the universe collapses rather quickly after it stops 
expanding. 

FigG] shows the constraints on the linear model from current observational data, assuming 
uniform priors on (p[^, Vq, a). The first row of figures in FigQ shows the probability distribution 
functions of the set of independent parameters Vq, p^, and a. The second row of figures in 
FigH shows derived parameters Q^, h, and the time to collapse from today t^. The dashed 
lines denote results using SN data only, the dotted lines denote results using SN together with 
CMB data, the solid lines denote results using SN together with CMB and LSS data. 

Table 1 lists the mean and the 68% and 95% confidence ranges of Qm, h, a, and tc from 
Figffl 
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157 SNe [Riess sample gold set] (dashed); SNe plus CMB (dotted); SNe plus CMB and 2dF (solid) 




[87rGV„/3](h/Ho)z [87TGp2/3](h/Ho)2 a[3/(877G)]'/2 




FIG. 1: The constraints on the hnear model parameters from current SN data only (Riess sample gold set, flux- 
averaged with Az ^ 0.05) (dashed hues), SN together with CMB data (dotted Hues), SN together with CMB 
and LSS data (solid lines). Uniform priors on {p'^„, Vb, a) are assumed. The first row shows the probability 
distribution functions of the set of independent parameters Vq, p^, and a The second row shows derived 
parameters f2„i, h, and the time to collapse from today tc- 



Table 1 



The constraints on the linear potential dark energy model from current data 





SN only'* 


SN plus CMB'' 


SN plus CMB and LSS'^ 




.25[.19,.30][.12,.35] 


.26 [.22, .30] [.18, .35] 


.26 [.22,.30][.19,.33] 




.651 [.642, .660] [.634, .668] 


.651 [.643, .660] [.635, .668] 


.652 [.644, .661] [.635, .668] 


" \ 8ttG) 


0.00 [-1.30, 1.30][-1. 60, 1.60] 


0.00 [-1.20, 1.20][-1. 60, 1.60] 


0.00 [-1.20, 1.20][-1.60, 1.60] 


tc/[Gyrr 


56.01 [35.75, — ] [19.18, —] 
[— , 108.18] [— , 2704.01] 


67.48 [42.95, — ] [23.61, — ] 
[— , 134.34] [— , 4095.18] 


65.66 [41.92, — ] [23.74, — ] 
[— , 130.52] [— , 3652.01] 


\minl -^dof 


25.98/23 


26.04 /24 


26.16/25 



" Riess sample gold set (157 SNe la), flux-averaged with Az = .05. 
^ CMB shift parameter R = ^llC^VizcMB) = 1-716 ± 0.062. 
The linear growth rate f{z2dF) = 0.51 ± 0.11. 

Statistical error only, not including the contribution from the much larger SN la absolute magnitude error of 
cr^"* ~ 0.05 Q. 

For tc, we list median, the 68% and 95% lower bounds and the 68% and 95% upper bounds. 
^ The number of degrees of freedom. 
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tc/to 



50 




0.25 0.5 0.75 1 1.25 1.5 



FIG. 2: The ratio of the collapse time from today, ic, and the age of the universe today, to, as a function of 
the linear potential parameter |a|. Fig^ 

Note that only the median, 68% and 95% confidence lower and upper limits of the collapse 
time from today tc are given, since a = is the cosmological constant model with tc —>■ oo (the 
mean of tc is not well defined for this reason). Computationally, we have to make a cutoff in 
for all models that are longer lived than the computation limit. The dependence of tc/to on a 
is shown in Fig|2l where to is the age of the universe today. 

Figini shows the effect of assuming different priors for a, using only SN data (Riess sample 
gold set, flux-averaged with Az = 0.05). The parameters are the same as in FigQ The solid, 
dashed, and dotted lines correspond to priors of p{a) oc 1, and a respectively. Table 2 

shows how assuming different priors for a changes the median, 68% and 95% lower and upper 
bounds on the collapse time from today tc- 
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157 SNe [Riess sample gold set], with prior p(a) «1 (solid); «l/a°* (dashed); «a (dotted); 




FIG. 3: The effect of assuming different priors for a, using only SN data (Riess sample gold set, flux-averaged 
with /S.Z = 0.05). The parameters are the same as in Fig^ The solid, dashed, and dotted lines correspond to 
priors oi p{a) oc 1, a~^^^, and a respectively. 

Table 2 

Effect of assuming different priors for a on the collapse time from today tc (in Gyrs) 



P{a) oc 


1 


^-0.5 


a 


median 


56.01 


96.92 


35.04 


68% lower bound 


35.75 


38.77 


27.41 


95% lower bound 


19.18 


19.38 


17.91 


68% upper bound 


108.18 


523.35 


48.75 


95% upper bound 


2704.01 


714260.81 


186.42 



IV. MCMC VERSUS FISHER MATRIX 



In this section we compare the joint confidence regions from our MCMC calculations to 
that of the popular Fisher matrix estimates. The Fisher matrix approach is an easy way to 
estimate the confidence regions in an n-dimensional parameter space pj, {i = 1, 2, ra), and it 
is especially powerful in the case when the probability distributions are close to the Gaussian 
distribution. The Fisher matrix is the inverse of the covariance matrix in this n-dimensional 
parameter space, see [46l |: 

F = C-\ 
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FIG. 4: The probability distribution function for the set of parameters \a\ and Qm with uniform prior on a 
(left figure) and with the prior proportional to a^'^'^ (right figure). 

It is easy to calculate the covariance matrix C using the MCMC data and marginalizing over 
additional parameters. We will consider the distribution and confidence regions for a and 
flm- In our case, the distribution of the parameter a (slope of the potential) is significantly 
non-Gaussian (see Fig^. There are a few facts that we have to take into account. First, 
the symmetry of the picture is dictated by the fact that the —a and a slopes in Eq.(l) are 
indistinguishable. It also leads to a paradoxical statement that Qm and a parameters become 
uncorrelated (the corresponding off-diagonal elements of the Fisher matrix are 0). Second, the 
actual parameter that we are trying to measure, and, hopefully, distinguish from (cosmological 
constant case) is actually This means that in our estimations for joint confidence regions 
we have to consider absolute value |a| rather than a. In Fig|3]the two-dimensional probability 
distribution functions for the set of parameters |a| and are presented for two different cases 
- SN plus CMB and LSS data with uniform priors on (p^, Vq, a) (right plot) and SN data with 
p{a) oc (left plot). In both cases the probability distributions are non-Gaussian. 

The confidence regions form the Fisher matrix and MCMC method are shown in FiglSl for 



FIG. 5: The confidence region for the set of parameters |a| and ilm- The purple contours are MCMC 68% and 
95% regions and red solid (dashed) ellipsoids correspond to 68% (95%) Fisher matrix estimations. 



0.46 




O.a 



a 
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SN with CMB and LSS data and uniform priors on {pl^, Vq, a) (left plot in FiglH). The MCMC 
confidence regions are presented by two uneven contours for 1-a (68%) and 2-ct (95%) values. 
Two Fisher ellipses (red solid ellipse corresponds to 68% confidence region and red dashed 
ellipse corresponds to 95% region) are shown as well. It is clear that due to the non-Gaussian 
nature of the distribution the Fisher ellipses could be used only for very rough estimations of 
the confidence regions. In fact, according to Figl^l a = (the cosmo logical constant case) 
would be excluded at 1-a level for the Fisher matrix approach, that is clearly not the case from 
the MCMC results. The 2-cr regions for both cases give similar results for a = and 10% 
discrepancy for the upper bound on 



V. PREDICTIONS FROM MODEL-INDEPENDENT BOUNDS ON THE PAST 
EVOLUTION 



We now derive complementary bounds on the cosmic doomsday by finding linear poten- 
tial models that fit completely under the 68% and 95% confidence contours of the model- 
independent reconstruction of dark energy density px{.z) from Ref . [lo| . 

The idea is to map out the boundaries of the past behavior of px{z) in a model- independent 
way based on SN la and other observations, in this case by the means of splines. Splines have 
the advantage over polynomials, for instance, that they do not have the tendency to wiggle in 
between measurement points. While the results will depend somewhat on the type of splines 
used, i.e. on the number of parameters (the more parameters one uses, the more information 
is extracted from the data, but the weaker the constraints will become), it is important to 
realize how this approach differs from earlier, conventional, model-dependent studies. The 
splines approach is able to model any arbitrary function, not making any assumption about 
the functional form of px{z) or wx{z)., contrary to the practice prevailing in present literature 
of using two parameter fits.^ 

As a starting point for our study, we take the model-independent boundaries on the past 
evolution of px{z) / px{Q) obtained in [l^. Next, we solve the differential equations for the 
evolution history of the universe with the linear potential for all possible initial conditions for 
the scalar field consistent with the 95% and 68% confidence regions of ^2^, and fit their cor- 
responding past px{z)/ px{0) history into the above envelope. The model parameters yielding 
the shortest lifetime, but still not being in violation anywhere with this envelope, determine the 
minimal lifetime obtainable by this method. These limiting cases, together with the boundaries, 
are graphically represented in FigEl 

Our approach is different from fl^ in that we recognize that an accurate prediction of the 
future can only be obtained by a model-specific extrapolation. One can study the past model- 
independently, yet whatever the conclusion from it is, the future will always depend on the 
specific underlying physical model. Even in the absence of knowledge of the exact physics — as 
it is the case presently with dark energy — one is forced to assume a specific physical model to 
obtain a prediction on the future lifetime of the universe. In our case, the model is the scalar 
field linear potential. 



3 Refs. [ill and |42| show the perils of such practice. 
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Px (z) /px (0) 




FIG. 6: The model-independent constraints on the past of px{z) / px{^) and how they constrain the hnear 
model. The bounds are represented by thick lines, while the corresponding linear model is shown by thin lines 
with same color and same texture (same dashing). Green denotes the 3-parameter spline bound and its linear 
model, while red shows the 4-parameter splines. The dashed lines depict the 68% confidence level, while the 
solid lines stand for 95%. Since the 3-parameter constraint contours always lie inside the 4-parameter spline at 
a given confidence level, one can argue that the 3-spline investigation is sufficient for this analysis, yielding a 
minimal future lifetime of 20.6 Gyr at 95% confidence level. 

We use the Hubble time H^^ as the time scale in our calculations. To translate the lifetime 
into Gyr, we need the value of the Hubble constant Hq. The Hubble constant, as obtained from 
the model-independent analysis in has a mean value of Hq = 63.7km/s/Mpc, corresponding 
to a Hubble time H^^ = 15.36Gyr, for the 3-parameter splines, while Hq = 63.5km/s/Mpc 
{Hq^ = 15.41Gyr) for the 4-parameter spline.^ 

Even the 95% error bars of the distribution for h do not cause an uncertainty of more than 
5.7%, corresponding to 0.87Gyr, thus we will neglect this uncertainty in our lifetime calculation 
in this section. 

In this manner, we obtain for the future lifetime of the universe: 



Table 3 

The lower bound on future lifetime of the universe 



Spline Type 


Confidence Level 


Future Lifetime 


4-parameter 


95% lower bound 
68% lower bound 


18.1 Gyr 

22.2 Gyr 


3-parameter 


95% lower bound 
68% lower bound 


20.6 Gyr 
25.3 Gyr 



These results turn out to be a little weaker than our model-specific analysis above. This 



^ These Hubble constant values are fully consistent with other direct measurements. SN la data typically give 
a Hq value consistent with Hq = 60 ± 10kms~^Mpc~^ 0. The Hubble Key Project (based on Cepheid 
distances) gives Hq = 72 ± 8kms^^Mpc~^ 0|. Note that CMB data do not give a direct measurement of 
Ho- 



is 



is not unexpected; it will be very difficult, if not impossible, to close the gap and raise the 
model-independent analysis to the same accuracy level as the model-dependent one. After all, 
no assumption about the model in the past is made. 

The model-independent analysis also depends on the number of free parameters used. The 
difference in results between the 3-parameter and the 4-parameter spline is roughly 10%. As 
usual, stronger constraints are achieved by fewer parameters at the cost of less flexibility. Since 
the 3-parameter spline contour lies completely within the 4-parameter one, one can argue that 
for the studied sample of supernovae using 3-parameter splines is sufficient. 

We conclude that the results obtained by the model-independent study of the past are some- 
what less stringent than our results obtained in Sec. 11111 Indeed, constraining the parameters 
of a given model directly is the most precise way to proceed. However, it is quite convenient 
to have a model-independent analysis of the past, and then apply it to any particular model 
for a quick estimate of the lifetime of the universe. This method provides a useful shortcut, 
which lets one to obtain approximate results without having to employ the whole machinery of 
MCMC that we have used to obtain our results in Sec. IIIII 



VI. DISCUSSION AND SUMMARY 

We have studied a representative dark energy model with a linear potential, V{(j)) = Vo(l + 
q;0), assuming a fiat universe. This model has the interesting property of dooming the universe 
to a Big Crunch for a ^ 0. It is the simplest doomsday model, in which the universe collapses 
rather quickly after it stops expanding. 

Since this model has only one free parameter, aVo, it is well constrained even by SN la data 
alone (Riess sample), see FigHJ 

We studied this model using current SN la (Riess sample), CMB (WMAP, CBI, and 
ACBAR), and LSS (2dF) data. We have used fiux-averaging to minimize the effect of weak 
lensing '14, 3f|, and included CMB and galaxy survey data in a consistent manner We 
found that in the context of the model with V{(f)) = Vo(l + acp), observations imply that the 
cosmic doomsday is more than about 42 billion years from today at 68% confidence, and more 
than about 24 billion years from today at at the 95% confidence level. 

These results represent a stronger bound on the doomsday time than the estimate of tc ^ H 
Gyr based on the previously available observational data |j30i| . When the data from SNAP and 
Planck satellites become available, it will be possible to further strengthen our bound up to 
tc > 40 Gyr at 95% confidence [30|. Note that these timescales have been derived assuming a 
uniform prior on a (the parameter of the linear model). 

When one evaluates the cosmological consequences of a given theory, one should specify not 
only its Lagrangian and the cosmological initial conditions, but also the theoretical expectations 
(priors) for the values of the parameters in the model. Since the change of the prior on the 
parameters is equivalent to changing the underlying physical model, there is nothing surprising 
in the dependence of the final results on the choice of the prior. For example, if one assumes, 
on the basis of some physical theory, that the cosmological constant is given by |A| = e'^, and 
all values of 7 = ln|A| are equally probable, one would come to a conclusion that the prior 
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probability for A is strongly peaked at A = 0. This would completely change not only the most 
probable value of dark energy density obtained from supernova experiments by using MCMC 
method, but also the famous anthropic bounds on A. However, in the absence of convincing 
theoretical arguments in favor of the uniform prior for In A, one usually assumes uniform or 
nearly uniform prior for A. 

Similarly, the main emphasis of our work was on the study of the linear model with the 
uniform prior for a. To make our results more general and robust, we have also considered 
other physically motivated priors on a as well, based on Changing the prior on a can 

change the bounds on the cosmic doomsday. We found that taking different priors, ~ a 
and ~ q;~°'^, altered the 95% confidence lower bound on obtained for = const 
only by about 10%. If one assume the uniforrn prior for In one would conclude that the 
universe will live for an exponentially long time (iol . 2^ 1 . However, if, just as in the case with the 



cosmological constant, one leaves theoretical speculations aside, one may want to concentrate 
on the results of our investigation with the uniform prior for a. 

Interestingly, we were able to obtain not only the lower bound on the doomsday time, but 
the upper bound as well. For the uniform prior on a the full bound on the doomsday time 
is 24 Gyr < < 3652 Gyr at the 95% confidence level (SN la plus CMB and LSS data, see 
Table 1). In other words, if one believes that the dark energy is described by the simplest 
linear model, and the parameter a can take any value with equal probability, then our universe 
most probably will not collapse earlier than in 24 billion years, but it will most probably live 
no longer than 3.65 x 10^ billion years. This doomsday prediction is mo del- dependent: it is 
a consequence of our assumption that dark energy is described by the linear model (1), and 
that all values of the parameter a are equally probable. On the other hand, as we already 
mentioned, this model is rather generic; it is the simplest representative of a broad class of the 
dark energy models where the effective cosmological constant can change continuously and the 
cosmological constant problem can be solved using anthropic considerations BHIillil. To 



improve our constraints on the cosmic doomsday time in this class of models one would need 
further cosmological observations and a deeper understanding of the theory of dark energy. 
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